Spin relaxation in an InAs quantum dot in the presence of terahertz driving fields 
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The spin relaxation in a ID InAs quantum dot with the Rashba spin-orbit coupling under driving 
THz magnetic fields is investigated by developing a kinetic equation with the help of the Floquet- 
Markov theory, which is generalized to the system with the spin-orbit coupling, to include both the 
strong driving field and the electron-phonon scattering. The spin relaxation time can be effectively 
prolonged or shortened by the terahertz magnetic field depending on the frequency and strength of 
the terahertz magnetic field. The effect can be understood as the sideband-modulated spin-phonon 
interaction. This offers an additional way to manipulate the spin relaxation time. 
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PACS numbers: 71.70.Ej, 72.25.Rb, 78.90.+t, 78.67.De 
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I. INTRODUCTION 

One of the goals of semiconductor spintronics 1,2 is to 
realize quantum information processing based on elec- 
tron/hole spins. Coherent oscillations of spin state driven 
by an AC magnetic/electric field, which is the key of such 
a goal, have been broadly studied. 3,4,5,6,7,8,9,10 ' 11 Exper- 
iments have demonstrated successfully the coherent os- 
cillations of electron spin by optical Stark effect in quan- 
tum wells on femtosecond time scale and by gigahertz- 
gate- voltage-controlled g-tensor modulation. 3,9 Recently, 
driven coherent oscillations of single spin in quantum 
dots (QDs) at hundreds of MHz has also been realized. 11 
Theoretically, Rashba and Efros showed perturbatively 
that electron spin can be manipulated by a weak in- 
plane time-dependent electric field via the spin-orbit cou- 
pling (SOC), such as the Rashba 12 and the Dresselhaus 13 
couplings. 4 They called it the electric-dipole spin res- 
onance (EDSR). Similar schemes have also been pro- 
posed in QDs. 7,8 Cheng and Wu have discussed non- 
perturbatively the effect of an intense terahertz (THz) 
electric field on two-dimensional electron gas (2DEG) 
with the Rashba SOC where the spin splitting is of the 
order of THz. 5 They showed that a THz electric field can 
strongly affect the density of states of the electron system 
and induce a THz spin oscillations. Similar effects have 
also been studied in QD by Jiang et al. 6 However, up 
till now there is no study on the spin dissipation effect, 
i.e., spin dephasing/relaxation in the AC-field driven sys- 
tems, especially from a fully microscopic approach. Re- 
cently Duckheim and Loss have studied the EDSR in 
the presence of disorder in a 2DEG, where the dissipa- 
tion effect is introduced by a relaxation time. 10 It has 
been shown in the system without the SOC that the dis- 
sipation under driving field can be very different from 
the one without the driving field. 14 ' 15 A full microscopic 
kinetic equation under driving field can be developed 
with the help of Floquet-Markov theory. 16 The Floquet- 
Markov theory combines the Floquet theory which can 



solve the time-dependent (periodic) Schodinger equation 
non-perturbatively with the Born-Markov approximation 
which is frequently used in deriving the kinetic equation 
with dissipations. 

In the present paper, we extend the Floquet-Markov 
approach to the system with the SOC. The change in 
density of states by the THz laser field 5 ' 6 implies a mod- 
ification of the spin related scattering and then the spin 
relaxation time. We study the spin relaxation in an InAs 
QD by developing a Floquet-Markov type kinetic equa- 
tion with the electron-acoustic-phonon scattering. The 
spin relaxation time is obtained by numerically solving 
the kinetic equation. We develop the model and for- 
mulism in Sec. II and present the numerical results in 
Sec. III. We conclude in Sec. IV. 



II. MODEL AND FORMALISM 

We construct our theory in a ID QD in which electron 
is strongly confined in the x-y directions and relatively 
weaker along the z-axis. This kind of QDs have been 
realized experimentally and are attracting more and more 
interests due to the good controllability of the size, shape, 
position, and electronic structures. 17,18 A static magnetic 
field and a THz electric/magnetic field, which can be 
provided by the free electron laser, 19 are applied along 
the z-axis. The total Hamiltonian is then given by 
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t 



H to t — H, 

t 



Hp ~\~ H e p , 
and H ep 



a 



(1) 

and the 



P ^qi? HV~qi)'~HV —™ — ejj ^— <qrj 

qj) )exp(iq ■ r) representing the phonon 
electron-phonon interaction Hamiltonian respectively. 
The electron Hamiltonian in the Coulomb gauge reads 20 
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+ V c (r)+H so (P)+H z 



(2) 



in which P = -ifiV + e/cA(t) with A(f) = 
— cEi cos(Oi)/0 + i[B + Bi cos(fli)] x r denoting the 
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vector potential induced by the THz electric/magnetic 
field and the static magnetic field, to* denotes the 
electron effective mass. V c (r) is the confining poten- 
tial of the QD which is taken to be square well with 
infinite well depth in each direction. 18 H so (P) is the 
SOC term, which consists of the Rashba term due to 
the structure inversion asymmetry 12 and the Dressel- 
haus term due to the bulk inversion asymmetry. 13 In 
InAs the Rashba SOC is the dominant contribution. 
Hz = 5<7/Ub(B + Bicos(fit)) • cr represents the Zee- 
man term. When the transverse confinement is strong 
enough, only the lowest subband is considered. Then the 
Hamiltonian of the electron can be reduced to an effec- 
tive Hamiltonian H eff = -^-r + V(z) + H*(P Z ) + H z , 



where P? 



Pz 



eE\ cos(fii)/fi and p z 
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Rashba term is Hf a — ~-Jf-(a x E y - a y E x )P z with E x 
and E y being the electric fields along the x- and in- 
directions which are due to the structure inversion asym- 
metry and can be controlled by the gate voltage. 21 For 
simplicity we choose E y = in our investigation and 
H so = -e-fE x a y P z = -f-<J y P z . This choice will not 
change the results of the calculation as one can always 
take a unitary transformation in the spin space to trans- 
fer any other configuration into our simplified one. The 
effective Hamiltonian is then written as 
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H eff (t)=H + H-(t)+H 2 (t) 



V(z) + -J^VyPz + ^9^bB <J z 
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Hi(t) = [\ 9 HBB-a z - + ^a v )\ cos(ffi) , (5) 



H 2 (t) = 



e 2 El cos 2 ■(tit) 
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Hence H eff (t+T ac ) = H eff (t) and T ac = ^ . We observe 
that H 2 (t) is only a function of time, and does not contain 
any other physical variables of electron. Thus it only 
induces a universal phase and has no contribution to the 
kinetics of the system. The corresponding Schrodinger 
equation 



iftjU(t) = H eff (t)*(t) 



(7) 



can be solved via the Floquct-Fourier approach developed 
by Shirley 22 and applied lately in systems with the SOC 
by Cheng and Wu. 5 The solution is 
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in which {(f> a (z)} is a complete set of the wave func- 
tions, chosen here to be the eigen- functions of a infinite- 
depth-square-well potential V(z). 5 ' 6 {e\} and {F* a } are 
the quasi-energies and the eigen-vectors of the following 
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where (r, t\an) = (j) a (r)e mVJ 
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H F = H eff (t) 



and {an\H r \l3m) = H^ m + mnS a0 S nm with H n 



tjt j 1 dte~ mSlt H e ff(t) representing the n-th Fourier com- 
ponent of the effective Hamiltonian. Due to the period- 
icity of Hjr, the eigenvalues are also periodic and can be 
written as s\j = 6\.q + Itt where e^.o is the eigenvalue in 
the region (— fi/2]. It is noted that s\j and £a,o cor- 
respond to the same physical solution to the Schrodinger 
equation. In the following we denote £a,o with ex for 
simplicity. The eigenvectors of the eigen-equations sat- 
isfy the orthogonal and complete relations 
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From these relations one obtains the orthogonal and com- 
plete relations of the Floquet wavefunctions 14 ' 22 

J2n(z,mx(z',t)=s(z- Z % (12) 

A 

/oo 
dz*l(z,t)* x ,(z,t) = 5x,x . (13) 
-oo 

The wavefunction [Eq. (8)] includes two significant ef- 
fects, one is the sideband effect 23 and the other is the AC 
Stark effect. 24 The former refers to the many frequencies 
s\ — nil in the wavefunction and the later represents the 
field-induced change of E\. 

The Floquet wavefunctions which contain all the dy- 
namic properties of the electron system without the 
electron-phonon coupling, give an optimal base to solve 
the equation of motion of the reduced density matrix of 
the electron system. The Floquet-Markov method which 
combines the Floquet solution of the electron Hamilto- 
nian and the Born-Markov approach to solve the equa- 
tion of motion under strong AC driving field was de- 
veloped by Kohler et al. in the absence of the SOC. 16 
Generally this method works well when the AC driven 
electron system is in the dynamic stable regime and the 
system interacts weakly with a Markovian reservoir with 
a damping rate much less than any eigen-frequency of 
the system. The latter requirement can be satisfied for 
almost every case in the spin decoherence problem, as 
spins are generally expected to have a very long coher- 
ence time. 25 In addition, studies have shown that, as well 
as keeping good quantitative results, the Floquet-Markov 
method has the advantage of being easy to handle nu- 
merically compared with the rather complicated path- 
integral approach. 16 Thus this method is very useful in 
the study of relaxation/dephasing in nano-structures. In 
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the present paper, we apply this method to the systems 
with the SOC to study the spin relaxation in QD due to 
the electron-acoustic-phonon scattering under the THz 
driving field. 

With the standard Feynman- Vernon initial condition 
p(to) = p e (to) <S> Peg, where p is the density matrix of the 

I 



electron and phonon system; p e is the density matrix of 
the electron system and p v eq represents the density matrix 
of the equilibrium phonon reservoir, and within the Born- 
Markov approximation, the reduced density matrix of the 
electron system satisfies the following equation: 



d i i r°° 

- p e = --[H e (t),p e ]- ¥ j dTTT p {{H ep ,[H ep (t-T,t),p e ®pr eq }}} (14) 

with Tr p standing for the trace over the phonon degree of freedom. H ep (t — r,t) = 11^(1 — T,t)H ep U~o(t — r,t), in 
which U a (t — T,t) = V t cxp[— jr f* T dt'(H e (t') + H p )] with V t denoting the time-ordering operator. Next we express 
this integral-differential equation of operators in a complete base of Floquet wavefunctions denoted by {\X(t))}. Using 
the complete relation J2\ | A(t)) (A(t) j = l, 14 ' 22 after some simple algebra, one has 

f) 1 poo 

A3 A4 

1 f'°° 

= 'ft ^EE(^A 3 ^Pi 4 A 2 -^A7A 3 Pi 3 A 4 ^A q 4 A 2 )((^W^-q,))+^. (15) 

A3A4 m 

where X^ iA2 = <A X (*)| exp(iq.r)|A 2 (t)), A m (t) = A^ta^e'^'+a^e-™') and X^ = (Ai(t)|£7 et (t-r, t) cxp(iq- 

t)Uq (t - r, t)\X 2 (t)) with U§(t — r,i) = V t exp[— ^ J t ~ T dt' H e (t')]. ((■ ■ ■ )) in Eq. (15) represents the statistical average 
over the phonon equilibrium distribution. By substituting the Floquet wave functions into X, one obtains AT q A = 

Ek^ iAxiX2kt Kx 2 k where A Al A 2 fc - (e Xl - e X2 )/h+kQ and X^ = J2 n J2 a p F na* F n+k ^(«|exp(iq ■ r)\(3) = 
^A 2 A*i-fc- As ^ ne Floquet wavefunctions are the solutions to Eq. (7), one has 

*lx 2 = (Ai(t)|C/ et (i-r,i)exp( iq .r)[/ e (t-T,t)|A 2 (i)) - (Ai(t - r)| cxp(zq • r)\X 2 (t - r)) 

= Y, eiAxiX2k(t ~ T)x lx 2 k- (16) 
k 

Therefore, Eq. (15) reads 

A 3 A 4 feife 2 q>) 

- X q X q * n\ \ P J ( A A 4 A 2 fc 1 -AA 3 A 1 fc 2 )t ( 7 (\ )\ + Hr (171 

vV A 4 A 2 fei" A -A 3 Aife 2 ^3A 4 e ^q>7^A 3 Aife 2 )} T n.C. , [II) 

with C qr) (A) = n(djq ri )S(A + w qr) ) + (n(u; qr) ) + 1)<5(A — w qj? ). Here n(w q?? ) is the Bose distribution function. The 
summations over k\ and fc 2 range from — oo to oo. The terms with C q?; (AA 3 Aifc 2 ) describe the /c 2 -photon-assisted 
scattering. These equations are still time-dependent. With the rotating wave approximation (RWA), one can sweep 
out the time-dependent terms which oscillate much faster than the damping rate of the density matrix 16 and Eq. (17) 
is further simplified into 

^PA X A 2 = - E AA i A 2A 3 A 4 pi 3 A 4 (18) 
A 3 A 4 

with 

A Al A 2 A 3 A 4 = {^2" E E^I M q r '| 2 E X A 1 A 5 fe 1 A "A 1 V5fe 2 5 fei,fe2^A 1 ,A 3 ^A 2 ,A 4 C' q , ) (AA 1 A5fc 2 ) 
fei,fe 2 <V1 A 5 

- X A 4 A 2 fc 1 X A 3 A 1 fe 2 ^A 4 -eA 2 -£A 3 +eA 1 ,(fc2-fci)^ C 'qr;( A A 3 A 1 fe 2 )]} + {Ai A 2 ,A 3 <-» A 4 j (19) 

I 

being a time-independent tensor in the case without de- generacy. {Ai <-> A 2 , A3 <-> A4}* in the above equation 
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stands for the same terms as in the previous {} but in- 
terchanging Xi and A 2 , A 3 and A4 and taking a com- 
plex conjugate. In the zero driving field limit, the Flo- 
quet wavefunction reduces from a many-frequency one 
to a single-frequency one. Thus, the only nonzero con- 
tribution in the summation X% ^ = e lAxiX2kt X^ iX2k 
is the term with k = k 0} with hAx 1 x 2 k = E\i — E\ 2 
being the energy difference between states |Ai) and 
exj+ejj.tfe-fcijn can be simplified into 
{5x 1 ,X 3 Sx 2 ,X 4 + S\ 1 ,x 2 6\ 3 ,\ i ). Therefore Eq. (18) can be 
simplified into 

§- t PU 2 = {-£ E i^i 2 i(Ai| e iqr iA 3 )i 2 p e Al A 2 

x [n((j qri )6(Ex 1 - E Xa + u qj) ) 
+(n(w q „) + l)S(E Xl - E x , - w q ,)]} 

+ {Ai^A 2 }* (20) 
for Ai 7^ A 2 , where E\ is the energy of the state |A) and 

q>j a 3 

X [n(w qj? )5(£;Ai - #A 3 + Wqr;) 

+ (S(w q ,) + 1)<5(-Ea! - ^a 3 - w qr) ) ] 
-^|(A 1 |e^|A 2 )|VA 2 A 2 +o; qr) ) 

A 2 

+(n(w q „) + l)S(Ex 2 - E Xl - WqiJ ) ] } . (21) 

These equations are consistent with the kinetic spin 
Bloch equations. 26 

Equation (18) can be rewritten in the matrix form as 
■§lP e — — A.p e , which is a standard first order differential 
equation. It can be solved through the eigenvalues and 
eigenvectors of the matrix A. Thus, for any observable 
O, 

0(t) = Tr(Op e ) 

= ]T (A 2 (i)|6|A 1 (i))P(A 1 A 2 )(A3A 4 ) 

xe " r(Wp i)W)^(°) ( 22 ) 

with r = P _1 AP being a diagonal matrix and P be- 
ing the transformation matrix. By solving Eqs. (8), (9), 
(19) and (22) numerically with an initial spin polariza- 
tion ^(O), one obtains the time evolution of S z . 

III. NUMERICAL RESULTS 

We consider an isolated ID In As QD, where the low- 
lying states can be approximated by eigenstates in an 
infinite- well-depth potential: 18 V c (r) = if < x < L a , 
< y < Lb and < z < L c and V c (r) = 00 elsewhere. 



We choose L a = L& = 20 nm and L c = 70 nm in the 
calculation. The separation between the first and the 
second subbands is about 9 meV (15 THz) along the 
^-direction, and 120 meV along the x (or y)-direction. 
By averaging over the lowest states in the x and y di- 
rections, one can turn the problem into an effective ID 
problem. We apply a static magnetic field B of 0.5/0.7/1 
T along the ,2-axis which corresponds to a Zeeman split- 
ting of about 0.4/0.56/0.8 meV (0.62/0.87/1.25 THz). 
In this energy range the electron-acoustic-phonon scat- 
tering is dominated by the deformation potential cou- 
pling in InAs. The corresponding scattering matrix reads 
Mqsi = S-y/ hq/2Dv s i, where H is the deformation poten- 
tial, D denotes the volume density and v s i stands for 
the longitudinal sound velocity. All the parameters used 
in the calculation are listed in Table I. 27 We take the 
cut-off frequency of the phonon reservoir to be the De- 
bye frequency ujd- The temperature is taken to be 200 
mK, corresponding to an energy of 0.016 meV (0.026 
THz), which is quite smaller than the other energies 
of the system (especially the Zeeman splitting energy), 
i.e., we study the spin relaxation in low temperature 
regime where the phonon-absorption processes are ener- 
getically unfavorable. The Rashba parameter is taken to 
be cur = 3.0 x 10~ 9 cV-cm. 21 By including all the scat- 
tering processes between the Floquet states due to the 
electron-phonon scattering, one can calculate the scat- 
tering matrix A [Eq. (19)]. With a preparation of occu- 
pying the first excited Zeeman state of H [see Eq.(4)] 
as the initial state, one can obtain the time-evolution of 
S z . By taking the envelope of S z and subtracting the 
cquillibrium spin polarization, we define T\ as the time 
needed for decay of the spin polarization by a factor of 
1/e. 



D 


5.9 x 10 3 kg/m 3 


9 


-14.7 


Vsl 


4.28 x 10 3 m/s 


COD 


32.7 THz 




5.8 eV 


* 

m 


0.0239 m 



TABLE I: Parameters used in the calculation 

At the very low temperature we study the spin relax- 
ation is due to the spin-flip transition between the low- 
est Zeeman sublevels. Recently Fonseca-Romero et al. 
studied a model two-level system coupled to an Ohmic 
reservoir via a x . 2S They showed that the pseudo-spin re- 
laxation and dephasing can be greatly modified by the 
driving field when it is in the type of Aa z cos(ftt). Their 
results show that at low temperature when the frequency 
is below the cut-off frequency of the reservoir, the driv- 
ing field enhances the pseudo-spin relaxation, otherwise 
impedes it. However, in QDs, spin is coupled indirectly 
with the phonon bath via the SOC. The effective spec- 
tral density of the spin-phonon coupling is generally not 
Ohmic. 29 Remarkably, in QD this spectral density can be 
controlled by the QD geometry (size, growth-direction, 
etc.), magnetic field, gate- voltage and the strength and 
symmetry of the SOC. 29 ' 30 ' 31 ' 32,33 ' 34,35 ' 36 ' 37 
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FIG. 1: (Color online) Spin relaxation time T\ as a function 
of THz magnetic field strength B\ for (a) Bo = 0.5 T and 
fl = 0.7 THz (green chain curve); B = 0.7 T and Q. = 0.9 
THz (blue dashed curve); Bo = 1.0 T and Q = 1.3 THz 
(red solid curve), all in resonance, and (b) Bo = 0.7 T, SI = 
0.45 THz (green chain curve); 0.9 THz (blue dashed curve, in 
resonance); and 1.8 THz (red solid curve). The black dotted 
curve in (a) is the spin relaxation time calculated from the 
simplified model Eq. (28) for B = 1.0 T and fi = 1.3 THz. 
(c) The spin relaxation time as a function of Bo without the 
THz magnetic field (black dotted curve), and results from Eq. 
(28) including all sidebands (red solid curve); only the 0, and 
±l-th sideband (blue dashed curve), and only the ±1 and 
±2-th sideband (green chain curve). Note the scale of Bo is 
on the upper frame of the figure. 



In our study, we investigate the spin relaxation in an 
InAs ID QD in the presence of a THz driving mag- 
netic field (no THz electric field included) along the 
z-axis (whose frequency is below the cut-off frequency) 
from a full microscopic approach. By including enough 
base-states and using the exact numerical diagonaliza- 



tion method, 30 we can calculate the Floquet wavefunc- 
tions. We then substitute these wavefunctions into the 
kinetic equations Eq. (18) to obtain the dissipative dy- 
namics of spin in QD under the driving field. We find 
that differing form the Ohmic spin-phonon coupling, 28 
the spin relaxation time here can be either prolonged or 
shortened by the sideband effect, depending on the fre- 
quency and strength of both the THz magnetic field and 
the static magnetic field, with the latter providing the 
Zeeman splitting. 

We plot the spin relaxation time T\ as a function of the 
THz magnetic field B\ in Fig. 1(a) with the electric field 
E = 0. The green-chain/blue-dashed/red-solid curves 
correspond to B =0.5/0.7/1.0 T and SI =0.7/0.9/1.3 
THz respectively. The THz magnetic field is tuned to be 
in resonance with the Zeeman splitting. It is seen from 
the figure that the spin relaxation time T% increases with 
the strength of the THz magnetic field B\, with some 
modulations. It increases more rapidly for the case with 
Bo = 1.0 T than the other two cases: at B\ w 5.6 T, T\ is 
increased about 25 times of the value at B\ = 0. It is also 
noted that due to the different modulations, when B\ < 1 
T, the spin relaxation time T\ decreases rapidly /mildly 
for B = 0.5 T/0.7 T but increases rapidly for B = 1.0 
T. We also plot T x at E = 0, and B = 0.7 T as a 
function of B\ with different frequencies (0.45, 0.9, and 
1.8 THz) in Fig. 1(b). Again the spin relaxation time 
increases overall with B\ but with different modulations. 
Nevertheless as a result of the modulation, when B\ < 
0.7 T, Ti decreases for the case with $1=0.45 THz and 
0.9 THz (mildly) but increases for the case with $1=1.8 
THz. 

These numerical results can be understood from the 
following simplified analysis. At zero electric field E = 
0, the effective Hamiltonian [Eqs. (3-6)] in the matrix 
element form reads 

(n,a\H eff \n',a'} = [e° + ^[E° z + E\ cos(nt)]a}s n , n , 
x5 a ^, + E so (n,a;n',a')5 { _ 1)rl+ril+1=Q 5 cr+al= o, (23) 



with 

tor of a 
ergy. 



z\n) = W-^-sin( 1 ^), and 



the eigen vec- 

E® = n 2 h 2 Tr 2 /(2m* L 2 ) is the subband en- 
E Q Z = gusBo and E], = gusBi represent the spin 
splittings due to the static and the THz magnetic field. 
E so (n,a;n',-a) = aa R (n\^p z \n') = o L ffi^ ly The 
last term of Eq. (23), which is the energy due to the SOC, 
is nonzero only for states with different spin and different 
parity. For very low temperature and for not too strong 
THz magnetic field, we can restrict ourselves to consider 
only the lowest two states due to the Zeeman splitting 
described by the time-independent parts of the effective 
Hamiltonian Eq. (23). These two states, denoted as \t; a ) 
with a =|, | (+, -), read 



|&>«|i,*> + £ 



E so (2n, -o-;l,cr) 
Ei — E 2n + E z a 



\2n,-a), (24) 
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from the perturbation. The matrix elements of the ef- 
fective Hamiltonian in the space of these two states are 
(€<r\H e ff\&) w ±[E Z + gfj, B B 1 CQs(Qi)](T6 <Ti<T i approxi- 
mately by retaining only the dominant terms. E z is the 
energy difference of the lowest two levels of H which 
is approximately gfi B Bo for not too large Bq. For this 
Hamiltonian the Schodingcr equation can be integrated 
out directly. The Floquet wavefunctions are therefore 

**(*) = cxp{-^[ — + — sin(fii)]}^ 



2h 

.ag^i B B t 

= cM - l —2fr~ 

x exp(— im£lt)€„ , 



2Mi 

,ag^ B Bi 



>2^ J ^ 2 mi ' 



(25) 



in which J m is the m-th Bessel function of the first 
kind and £ CT represents the pseudo-spin wavefunction 
of the lowest two states [Eq. (24)]. The form of 
the wavefunction clearly indicates the sideband effect. 
Specifically, the probability of finding the electron in 
is J m ( gg 2^ Bl )l 2 which oscillates at a frequency 
ag^i B B /h + mfl. The frequency together with its cor- 
responding coefficient J m ( <79 2^f 1 ) is referred to as the 
m-th sideband from now on. In the following we show 
that this sideband effect can greatly affect the spin re- 
laxation. 

Substituting these wavefunctions into the kinetic equa- 
tion Eq. (17), one obtains the following scattering matrix 

Atttt = §f EI^I'E^tVI 2 ^^) 

qr; k 

= E^(^Tlfc) E J-mJm+kJ-m'Jm'+k(^^^Q ) 
k m,m' 



(26) 



Here the relation ^2 m J- m J m +k{x) = Jki^x) is applied 
and 

r(A nfe ) = |j^|M q7) | 2 |^ T |e^ r |^)| 2 [n(a; q7) ) 



xS(A nk + w q7) ) + (n(w q77 ) + l)5(A nfc - w q „)] . (27) 

Without the driving field, Eq. (26) reduces back to the 
well known form A-j-f-f-f = r(Afj,o)- It is noted that un- 
like the driving-field free case where only phonons with 
energy g/i B Bo contribute to the spin relaxation, due to 
the sideband effect caused by the driving field, from Eq. 
(26) one can see that phonons with energies Afj^ = 
g^isBo + kMl (k ^ 0) also contribute to the spin relax- 
ation. The spin relaxation rate is T-f 1 = 2Af jjj. Partic- 
ularly, at zero driving field, Tf 1 = 2r( 5 /x B £ )- The spin 
relaxation time as a function of the THz magnetic field 



B\ and the static magnetic field Bq, Ti(B\,Bq), reads 
Tx{Bx,Bq) » T 1 {0,B )/{j${g l i B B 1 /hn) 

+ £ )Jfe2( ^lr 1) r„ JSo B +khn) I (28) 

/c— — oo k—1 

in which T (gfi B B a + kMl) is the spin relaxation time 
at same external condition (QD geometry, static mag- 
netic field, temperature, etc.) but at different energy 
g[i B Bo + kMl. We approximately take To(g[i B Bo + kMl) 
as the spin relaxation time Ti(0,Bq), where Bq is de- 
termine by the condition that the corresponding Zeeman 
splitting for the lowest two states is gfi B B + kMl. This 
approximation is reasonable when g[i B Bo + kMl is much 
smaller than the energy difference between the first and 
second subbands, since the difference in the spin mixing 
for the lowest two states between the case with a static 
magnetic field Bq and the case with B is marginal. 

We first calculate 2i(0, -Bo) as a function of B and 
plot it in Fig. 1(c) as black dotted curve. Facilitated 
with this quantity, we further obtain Ti(Bi,B ) in Eq. 
(28) as a function of B x for Bq = IT and SI = 1.3 THz 
and plot it as black dotted curve in Fig. 1(a). It is seen 
that this approximate results agrees with the numerical 
result qualitatively. It is seen from the simplified model 
that the dominant contribution comes from the term of 
J o( ° 9> m Bl ) in thc denominator in Eq. (28) which is zero 
at B\ «2.4 and 5.6 T, corresponding to the first and 
the second peaks in Fig. 1(a). In order to reveal which 
sideband contributes to the peaks, we plot T\{B\,B§) 
calculated from Eq. (28) with all the sidebands [red solid 
curve, same as the black dotted curve in Fig. 1(a)], with 
only the 0, ±l-th sidebands (blue dashed curve), and 
with only the ±1, ±2-th sidebands (green chain curve) 
in Fig. 1(c) versus B 1 when B =l T and fl = 1.3 THz. 
It is seen in the figure that the result with only thc 0, 
±l-th sideband agrees with the total approximation re- 
sult pretty well when B\ < 3.5 T which indicates that 
the first peak is mainly due to these sidebands. When 
Bi > 3.5 T, the result including only the ±1, and ±2- 
th sidebands is in reasonable agreement with the result 
with all thc sidebands. Therefore the second peak mainly 
comes from the ±1 and ±2-th sidebands. The oscillations 
in spin relaxation time are mainly due to oscillations of 
the sideband amplitude with the strength of the driving 
field, i.e., the sideband factor J^.(gn B Bi/Ml) in Eq. (28). 
The relaxation time To(gfi B Bo + kMl) can be larger or 
smaller than Ti(0, Bo) depending on the details of the 
SOC mediated spin-phonon coupling, 29 ' 30 ' 31 - 32 ' 33 - 34 ' 35 ' 36 
which determines the effect of the THz magnetic field on 
the spin relaxation time. 

At very low temperature, if g/j, B B + kMl is less 
than zero, thc relaxation process is prohibited be- 
cause there is no phonon to be absorbed. There- 
fore, Ti(0, B )/T (gfj, B B + kMl) is zero for gfi B B + 
kMl < 0. For the case of resonant driven system, 
only k > terms remain finite. Among these terms 
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only the terms with k = and 1 are important for 
small B\. J^guBBx/m) w 1 - 2Jf(gf J , B B 1 /Ml) if 
gji B B\/Ml < 1. Thus the denominator is approxi- 
matly 1 + ^(gfigB^m) [7i(0, B )/T {gn B B + Ml) - 2] . 
Therefore, T (g[i B Bo + Ml) becomes an important factor 
for T\(Bi, Bo). The three color curves in Fig. 1(a) with 
B = 1.0/0.7/0.5 T, corresponding to T 1 (B 1 ,B ) be- 
ing increased/insensitive/decreased with B\ when £?i < 
1.1 T. This is because T (gfi B B + Ml) w Ti (0,2.80) 
is larger than/approximately/smaller than Ti(0,B )/2. 
The same thing happens, if one changes the frequency 
of the THz magnetic field, as To(gn B Bo + M7) also de- 
pends on Q. The property of the frequency dependence 
in Fig. 1(b) when B\ < 0.7 T can be understood from 
the fact that T (gii B B + Ml) « Ti(0,B + hfl/g^s) 
is smaller/a little smaller/larger than Ti(0,B )/2 when 
O = 0.45/0.9/1.8 THz. This indicates that by prop- 
erly tuning the frequency of the THz magnetic field, 
we can change the effect of the THz magnetic field effi- 
ciently. For strong THz magnetic field, the the sideband 
factor J^(gfj, B B 1 /hil) is important only for terms with 
large k where T (gfi B B + kMl) w Ti(0, B + kMl/g^ B ) 
is larger than Ti(0,i?o) for large enough k. Moreover, 
J%(g/j, B Bi/hn) also decreases with B\. Therefore for 
strong enough THz magnetic field, the spin relaxation 
time is always larger than Ti(0,B o ) as indicated by the 
six colored curves in Figs. 1(a) and (b). 

IV. CONCLUSIONS 

In conclusion, we apply the Floquet-Markov theory to 
the spin kinetics in ID InAs QD to study the spin relax- 



ation in the presence of the THz driving field. Especially, 
we study that the spin relxation under a THz magnetic 
field which is parallel to a static magnetic field. We find 
that the spin relaxation time can be effectively manipu- 
lated by the driving field depending on its frequency and 
strength. This offers a new way to control the spin re- 
laxation. The effect is understood as the sideband effect 
modulates the indirect spin-phonon coupling. The effect 
of the driving field also depends on the properties of the 
QD, such as the QD geometry, the strength and symme- 
tries of the spin-orbit coupling, etc., which can be tuned 
by the gate-voltage, and the static magnetic field. The 
formulism developed here can be generalized to other sys- 
tems, such as the two-dimensional electron/hole gas with 
the SOC to study the spin relaxation. The corresponding 
kinetic equation can further be used to study the prob- 
lems such as the AC-field-induced spin polarization and 
the related spin transport. These are still under investi- 
gation and will be published elsewhere. 
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